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AIR FORCES AND MOMENTS ON TRIANGULAR AND RELATED 
WINGS WITH SUBSONIC LEADING EDGES OSCILLATING 
IN SUPERSONIC POTENTIAL FLOW 
By Charles E. Watkins 

SUMMARY 


This analysis treats the air forces and moments in supersonic 
potential flow on oscillating triangular wings and a series of swept- 
back and arrow wings with subsonic leading edges and supersonic trailing 
edges. For the wings undergoing sinusoidal torsional oscillations simul- 
taneously with vertical translations, the linearized velocity potential 
is derived in the form of a power series in terms of a frequency param- 
eter. This method can be useful for treatment of similar problems for 
other plan forms and for wings undergoing other sinusoidal motions. 

For triangular wings, as many terms of such a series expansion as may 
be desired can be determined; however, the terms after the first few 
become very cumbersome. 

Closed expressions that include the reduced frequency to the third 
power, an order which is sufficient for a large class of practical 
applications, are given for the velocity potential and for the components 
of chordwise section force and moment coefficients. 

These wings are found to exhibit the possibility of undamped tor- 
sional oscillations for certain ranges of Mach number and locations of 
the axis of rotation. The ranges of these parameters are delineated for 
triangular wings. 


INTRODUCTION 


This paper is concerned with the derivation of expressions for the 
velocity potential and associated forces and moments for oscillating 
triangular wings in supersonic flow. The boundary- value problem for 
the linearized velocity potential for an apex-forward triangular wing 
oscillating in a supersonic main stream may be classified, according to 
reference 1, as "purely supersonic" if the leading edges of the triangle 
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are outside the Mach cone emanating from the apex of the triangle or 
"mixed supersonic" if the leading edges are inside this Mach cone. 

In the purely supersonic case the principle of independence holds; 
that is, the flow on the upper surface of the wing is independent of the 
flow on the lower surface and vice versa. Garrick and Rubinow (refer- 
ence l) have shown that the boundary -value problem for the velocity 
potential in the purely supersonic case can be satisfied by simple dis- 
tributions of sources with local strength proportional to the local 
prescribed normal velocity of the wing. 

In the mixed supersonic case the principle of independence does not 
hold. Boundary -value problems for lift-producing wings in this case can 
be satisfied by distributions of doublets; the relation between doublet 
strength and normal velocity of the wing is, however, in general, not 
simple. The determination of this relation requires the solution of an 
integral equation that employs the potential of a time-dependent unit 
doublet as kernel and limits of integration that depend on Mach number 
and wing plan form. 

For treatment of problems that involve boundary conditions that are 
independent of time, such as constant angle of attack, constant rate of 
pitching, and so forth, the doublet potential, employed as kernel of the 
integral equation, is considered independent of time and in these cases 
the integral equations for triangular wings can be solved by a straight- 
forward process. 

For treatment of problems of oscillating wings, however, it is 
necessary to employ, as the kernel of the integral equations, a doublet 
potential that varies harmonically with time and in this case the solu- 
tion of the integral equation, generally, becomes very cumbersome. If 
the doublet potential or kernel is expanded in terms of the frequency 
of oscillation, however, use can be made of knowledge of solutions of 
integral equations for problems that are independent of time to obtain 
an expanded form of solution for a wing undergoing harmonic oscillations. 
Such a procedure was demonstrated in treatments of rectangular wings in 
references 2 and 3 * (The derivation in reference 3 is based on an 
erroneous argument, regarding certain terms in the normal velocity; 
nevertheless, the final expression given for the velocity potential is 
correct. ) 

The purpose of the present paper is to make use of the expanded 
form of the velocity potential to obtain the forces and moments, based 
on the first few terms of this potential, for a rigid triangular wing 
performing vertical and pitching sinusoidal oscillations in mixed super- 
sonic flow. Although as many terms of the expanded potential as may be 
desired can be obtained after the first few terms, the process becomes 
very cumbersome. The flow normal to the leading edge is subsonic but 
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the flow normal to the trailing edge is considered to be supersonic. 

This latter consideration implies that the potential derived for trian- 
gular plan forms may be used to calculate the aerodynamic forces and 
moments for other plan forms that may be formed with the triangular wing 
by cutting the trailing edges so that they lie ahead of the Mach cones 
emanating from their foremost points. 

Other approaches to the solution of the problem of oscillating 
triangular wings have been given by Robinson (reference 4), Haskind and 
Falkovich (reference 5), and by Stewartson (reference 6). In both 
references 4 and 5 formal solutions to the problem were obtained in terms 
of special systems of curvilinear coordinates. Robinson's solution was 
given in terms of a double summation of trilinear combinations of Bessel 
functions of the first kind with Lame' - functions of the first and second 
kinds. Similarly, the solution of Haskind and Falkovich was given in 
terms of summations of Bessel functions of the first kind combined with 
elliptic integrals of the first and second kinds. In both references 4 
and 5 the potentials were not reduced to useful forms for calculating 
forces and moments. 

In reference 6 Stewartson makes an interesting though specialized 
use of the Laplace transformation to develop a method whereby terms of 
the velocity potential for triangular plan forms, expanded as herein, 
can be obtained. Stewartson gives formulas that, except for errors 
presumably in printing, can be used to develop the potential to the 
second power of the frequency but he omits many details in his derivation. 


SYMBOLS 


0 disturbance-velocity potential 

x,y, z rectangular coordinates attached to wing moving in 

negative x-direction 

i ^ T| rectangular coordinates used to represent space location 

of doublets in xy-plane 

^ function defining mean ordinates of any chordwise section 

of wing such as y = y^ as shown in figure 1 

w(x,yi,t) vertical velocity at surface of wing along chordwise 

section at y = yp 


*0 


abscissa of axis of rotation of wing as shown in figure 1 
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t 

time 

h 

vertical displacement of axis of rotation 

*0 

amplitude of vertical displacement of axis of rotation, 
positive downward 

a 

angle of attack 


amplitude of angular displacement about axis of rotation, 
positive leading edge up 

• • 

h,a 

time derivatives of h and a, respectively 

V 

velocity of main stream 

c 

velocity of sound 

M 

free-stream Mach number (v/c) 


3 = V M 2 - 1 


0) 

frequency of oscillations 

- M 2 o> 

^ = — o 

V3 2 


k 

reduced frequency 

6 

half apex angle 


C = tan € 


a nm 

represents functions of ao, x, and M 

V F n 

functions used to denote doublet distribution functions 

d n 

constants associated with D n depending on PC 

A i,Xj 

constants depending on 3C 

°k 

constants depending on 3C and M 

2b 

root chord of wing 
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s 

0 

q,N,0,T 

P 

Ap 

P 

^1 


semispan of wing 

slope of ray passing through vertex of wing 

dummy variables 

density 

local pressure difference 

section force (total force at any spanwise station) 
components of section force coefficients 


M a section moment (total moment about x = xq at any 

spanwise station) 

components of section moment coefficients 


% 


total component of damping- moment coefficient 


ANALYSIS 

Boundary-Value Problem for the Velocity Potential 


Referred to a rectangular coordinate system moving forward at a 
uniform supersonic speed in the negative x-direction (see fig. l) the 
differential equation for the propagation of small disturbances that 
must be satisfied by the velocity potential is 


1 b 






+ 

^ Sz^ 


( 1 ) 


The main governing boundary condition to be satisfied by the velocity 
potential is that the flow be tangent to the surface of the wing, or 



— >0 


w(x,y,t) 


SZm SZjh 
Sx St 


( 2 ) 


where is the vertical displacement of any point of the wing. For 

the particular case of a wing independently performing small sinusoidal 
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pitching oscillations of amplitude Oq about some spanwise axis Xq 
and small sinusoidal vertical translations of amplitude the 

quantity in equation (2) is 

Zm = e tot |a 0 (x - Xq) + hj] = a(x - Xq) + h (3) 

(See fig. l(b) for sketch shoving instantaneous displacement of sec- 
tion y = y 1# ) For convenience, the frequency of oscillation of both 
pitching and translatory motion is denoted by a). Considering these 
motions to occur at separate frequencies would add no diff iculties to 
the derivation. 

Substituting the expression for (equation (3)) into equation (2) 

gives 

v(x,y,t) = Va + a(x - xq) + h (4) 


Equation (^-) implies that the velocity potential may be expressed as the 
sum of separate effects due to position and motion of the wing associated 
with individual terms of this equation, namely 


0 = 0a + Ax + ^h 


(5) 


Derivation of 0 

In order to obtain the analytical expression for the potential 0, 
it is necessary to derive only one of the subsidiary potentials appearing 
in equation (5), say 0 a . The other subsidiary potentials 0^ and 0^ 

can then be obtained from the derived expression for 0 a by simple 
comparison. 

In order to satisfy the boundary- value problem for 0 a , a convenient 

procedure is to start with the expanded form of the potential of a uni- 
form distribution of doublets. Then, for a given power of the frequency 
of oscillation this potential, as will be shown in the following analysis, 
can be modified so that, when integration is made over the appropriate 
region, the results satisfy the differential equation (l) to the given 
power of the frequency and satisfy the condition of tangential flow 
exactly. The type of doublet required is that with its axis normal to 
the plane of the wing. The potential of such a - doublet may be obtained 
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from the potential of a source, located in the plane of the wing, by 
partial differentiation with respect to the direction normal to the 
plane of the wing. Similarly, the potential of a distribution of the 
required type of doublets can be obtained from a distribution of sources 
located in the plane of the wing. 

The potential at (x,y,z) due to sources located at points (|,tj, 0) 
in region r (illustrated in fig. 2) of the xy-plane which satisfies 
the differential equation (l) may be written as 


Expanding the integrand of equation (6) into a power series in in, 
collecting terms with respect to |, and differentiating the resulting 
integral with respect to z gives the expanded form of the potential 
of a uniform distribution of doublets, namely 



( 6 ) 


where 


Mjo _ M^cjd 

cp 2 vp 2 


and 


r = ^(x - O 2 - P 2 (y - n ) 2 - P 2 z 2 



r 




at dq 


( 7 ) 
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where 


_ (M n O iok 


n! 


( 8 ) 


and 


\m-l /— \2m-2 


±i r^js\ a 

= (2» - 2) Am 1 111 


( 9 ) 


For convenience in the succeeding discussion and analysis equa- 
tion (7) may be written in the following form: 


= V" 


t n 


anl / / Y dl dT1 + 


£ fe]T ^ in 

(2m - 2)1 W ^ J r J 


(10) 


An interesting and significant property of equation (10) is that 
the coefficient of each power of I satisfies the differential equa- 
tion (l) and has the form of a source potential with strength propor- 
tional to This property may be shown by writing the coef- 
ficient of | n as follows: 


fnl 

R 


1 + 


m= 



2m-2 

fe") 

1 


nl / 0) _ 
cos — R 
M 


■( 


— \n _ 

n! ) 6 


-ico(x-£) 
e v ' cos 


1 ») 


( 11 ) 


2 


MCA TN 2^57 


9 


A more general solution to equation (l) may thus be obtained by intro- 
ducing properly chosen weight or distribution functions (denoted by 
£D n ( £> T l) ) into the coefficients of £ n in equation (10). Let this 
solution be denoted by 0^.1 then it can be written as 


00 



^n+1 

R 


d| dT] + 


oo 


I 

m=‘ 


(-D 11 - 1 

(2m - 2 ): 


CD 


2m-2 

a nl 



D n U,T)H n+1 R 2m 3 d| dt) 


( 12 ) 


Examination of equation (12) shows that, at the surface z = 0, the 
potential 0p is determined by the first integral expression but that 

both integral expressions may give rise to normal velocity. In succeeding 
steps in this analysis it is shown that the distribution functions D n (|,T]) 

in equation (12) can be determined so that the first integral expression 
taken alone will exactly satisfy the boundary condition of tangential 

/d0 \ 

flow for 0 a ; that is j = Va. Also, any additional normal veloc- 


ity that arises from the second integral expression can be canceled, to 
the required order, by consideration of additional doublet solutions to 
equation (l). The problem of satisfying the boundary-value problem for 
the velocity potential 0 a may thus be reduced to that of determining 

the appropriate distribution functions and additional solutions to 
equation ( l) . 


In order to show that the first integral expression in equation (12) 
can be made to satisfy the boundary condition for 0 a , the coeffi- 
cient Vaa^ appearing in this equation is first considered. If the 
analytical expression for the coefficient a^ (equation (8)) is multi- 
plied by x n and summed with respect to n, the result is identically 
Va. This result may be shown as follows: 

00 

Va 

n=0 


nl 


= Va 


00 

r 

n=0 


( icDX ) n -idix 


= Vae 


icjox -ia)x _ 


= Va 


n: 


( 13 ) 
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It is significant that this identity holds if only terms in o3 to any 
given power are considered. For example, retaining only terms including 
S3 to the third power gives 


Va 



icnx 


a^x^ ia>3x3 

“2“ + ~ 5~ 


+ x id) + 


S?x 




Va 


(14) 


Next consider the normal velocity at z = 0 associated with namely 


W 1 = 




tn+1 

Dn(l^) — ^ *) + 


V 2 - (-l) 111 ' 1 

^2(2m-2)!\M/ &nl 



D n (i,Ti)i n+;L R 2m ‘ 3 d| dn 

(15) 


Examination of this equation and equation (13) reveals that the first 
integral expression on the right of equation (15) yields Va exactly, 
provided the distribution functions are determined so that the following 
integral equation is satisfied: 


Va lim d 2 

n z >0 ^ 



^n+1 

R 


d| 


dq = Vax 11 


(16) 


The kernel of this integral equation has the form of a steady-state 
doublet potential. The problem of determining the distribution functions 
for this case is therefore analogous to determining distribution functions 
for certain steady-state problems. The distribution functions for steady- 
state problems, at least for those involving conical flow, can be deter- 
mined by a straightforward process, the main details of which are given 
in the appendix. In this appendix a method of solving equation (l6) for 
a triangular wing is derived and the distribution functions required to 
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derive the velocity potential for this wing to the third power of 
are given. It is to be noted that the method derived for solving equa- 
tion (l6) for a triangular wing may be generalized to apply to various 
plan forms and to problems of satisfying the boundary conditions for 
various velocity distributions. 

From this point on, the analysis is restricted to the derivation of 
terms of the expanded potential involving 00 to the third power. The 
method for deriving these first few terms is quite general and can be 
used to obtain as many additional terms of the expanded potential as may 
be desired. As previously pointed out, however, terms of the potential 
after the first few become very unwieldy. 

If the_appropriate distribution functions are known for terms 
involving co to the third power, equation (15) may be written as follows: 


Wj = Va + w 2 


(17) 


where 

Va lim d 2 a£ f f [7 _ . _ 0 ~| 

W 2 ‘ — z— >0 ^2^ J r J P " ^ D 0 R + ia ^ D 1 R J dl dq (l8) 


is the additional vertical velocity arising from the second integral 
expression in equation (12) involving d> to the third power. In order 
to maintain the boundary condition for 0 a , this additional velocity w 2 
must be canceled. As previously pointed out this canceling, to a 
required order, can be achieved by considering other doublet solutions 
to equation (l) . For this particular case consider a relation 0-^ 

similar to (equation (12)) having the following form: 


h 


Va a?_ 
n dz 2M 2 




Dn(l>n)l n+3 clS dl) + 


y - M )” -1 tif *- 2 

m=2 (2m - 2)!\M/ &nl 



D n ( |,ti) | n+ 3R 2m- 3 d | dT] 


(19) 
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The vertical velocity w^, to the third power of co, arising from this 
expression is 


W 1 


Va lim /a)\ 2 
2« z — >0 ^ z 2 \M/ 



(l - 1 cox)Dq 


ll 

R 


+ in>D 



d| drj 


( 20 ) 


In this equation the distribution functions Dq and Dp can be deter- 
mined, as discussed subsequently, by the method given in the appendix so 
that wp is identically equal in value but opposite in sign to w 2 . 

When these functions are determined, the boundary condition for 0 a is 
satisfied by wp + wp = Va, which implies that the potential 0 a to the 
third power of 53 is given by the sum 


fa = h + 01 ( 21 ) 

to this power of an 

Expressions for Dq and Dp are given together with expressions 
for other distribution functions Dq, Dp, D 2 , and D^ in the appendix 
In regard to the determination of Dq and Dp, as well as other orders 
for D n in more extended treatments, it may appear necessary, in order 

to formulate integral equations for these functions, to perform the 
generally unwieldy integrations of the type 


f f V *0 < 22 > 

r 

appearing in equation ( 15) • In. general, however, the information neces- 
sary for the determination of the functions D n can be obtained, as is 
done in the derivation of the functions D n in the appendix, by exam- 
ining the values of these integrals and their derivatives with regard to 

the parameter 6 = ^ at some particular value of 0 . 

Returning to equation (21) and introducing into this equation the 
expression for 0p (equation (12)) and the expression for (Dp 
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(equation (l9)) 7 each to the third power of gives for the 
potential 0 a to this power of ao the following result: 


f 


/yl 

S’ 

( 

'r J 

\ 


1 - fix - ^ X 2 + ^ x 3 \i . . ^4 xIR 

2 6 R \2M 2 2M 2 


| 2 Di 


^ . —2 ixo3 2\l ia>3 

iCD + - x I— - — — R 


2 R 2M 2 


- |3d 2 /(0 ‘ 


2 


iaj3 \l 

3 X /R 


t 4_ ia)3 1 3 - / co 2 ia)3 \l 4^ ia)3 1 

1 3 6 R 1 ^°\2M 2 " 2M 2 > 6 1 2M 2 R 


d| dq 


(23) 


Since it may be shown, as in reference 7* for example, that 

T I T / / D » |,1+1 « dT > * 

r 

equation (23) reduces, at z = 0, to 


0a - Va 


fl " IflK “ ~ x 2 + ^ x3)D Q (x,y) 


-2 iaj3x 2 

x ia) + (d x - 


■jl>l(x,y) - x2 (y - xjD 2 (x ;> y) - 


x 3 D 3 (x,y) + x 2 ^ - x)D 0 (x,y) + x3 ^ D^x^y) 

6 3 \2M 2 2M 2 / 2M 2 


( 24 ) 
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which, after the expressions for the distribution functions given in the 
appendix are substituted and the terms are regrouped, may be written in 
the following simple form: 


The quantities A q, Aj_, and a 3> and. a 4 i n "this equation 

are part of a group of quantities A-j_, Xj, cr k , which are defined in 

the appendix; these quantities are suitable for writing the potential 
and expressions subsequently derived for forces and moments in simple form. 
The quantities A^ and Xj are functions only of the product pc (ratio 

of tangent of the half apex angle of the triangle to the Mach angle) and 
are shown plotted in figures 3 and 4, respectively. The quantities cr k 

are functions of Mach number M and the product PC. These quantities 
may be evaluated for particular values of M from the plots of in 

figure 4. 

The quantities Aq and A^ are the same, as should be expected, 
as the parameters associated, respectively, with constant angle of attack 
and constant pitching of triangular wings of references 8 and 9* 

Expressions for the potentials 0^ and 0^ can be obtained by the 
method discussed for obtaining 0 a , or they can be obtained to the third 

order of cu by comparison and synthesis from equation (23). After 
simplification these expressions are 



(25) 



(26) 


K = h ]/c2 



^Aq - 2X-l1o)x 
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At PC = 1 or C = which is the condition at which the Mach lines 

from the apex of the triangle coincide with the leading edges of the 
triangle, equations (23)* (26), and (27) reduce, respectively, to 


K = ~\R - fV(i - 


icpx 7M^ + 
3 90M 2 


— CJD 2 X 2 


sW + 

90M 2 


2?^Ii<s3x3 

210M 2 


2M - 7 

210M 2 


io^P^y^x j 


(28) 



8ico x 2 + 2ico p2y2 _ ^(llM 2 + 7) x 3 + 
^5 ^5 315M 2 


a^(2M 2 + 7) 

315M 2 




io>x 7M 2 + 5 -2 2 2M 2 - 5 

3 90M 2 90M 2 



(29) 


K ■ j|V x2 - - t 1 - ^ x2 ' ^ 2 ) 


90M 


(30) 


For M = 1 or for values of the product PC such that P 2 C 2 « 1, 
equations (23)* (26 ), and (27) reduce, respectively, to 


0a 


Ycojc “X' 


2„2 



(31) 


= a-vt 2 ^ 2 - y 2 (x - Xq^ 


(32) 
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K = 



( 33 ) 


These expressions are the counterpart of the potential for vanishingly 
small triangular wings in steady flow (reference 10) and may thus he 
regarded as the potential of a small triangle oscillating in pitch and 
vertical translation in either subsonic or supersonic flow. 


As pointed out in the introduction, the velocity potential for the 
triangular wing can be used to calculate the aerodynamic forces and 
moments for other plan forms that can be formed from the triangular wing 
by cutting the trailing edges so that they lie ahead of the Mach cones 
emanating from their foremost points. Sketches of different plan forms 
thus obtained are shown in figure 5 • 

The force and moment coefficients desirable for most flutter cal- 
culations are those that yield the spanwise variation in these quantities 
or chordwise force and moment coefficients. These coefficients are 
obtained by integrating the pressure difference along any chord for the 
forces and the pressure difference multiplied by a moment arm for the 
moments. A convenient procedure in deriving these quantities is to 

introduce the reduced- frequency parameter ^ = k and to employ the 

the variables x, y , and xq in a new sense as nondimens ional quan- 
tities obtained by dividing the old variables by the maximum chord 2b 
of the wing. 

The pressure difference between the upper and lower surfaces of the 
wing is 


The force, positive downward, at any section of any of the plan 
forms shown in figure 5 may be expressed as 


Forces and Moments 



(3b) 
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P = 2b / Z!j> dx = -4pb 

y/c 


1 / v M + ^ 


y/c 


1 2b 


dx 


= -2pV J 0 


X 1 p x 1 . 

+ 2ik / 0 dx 

y/C U y/c 


(35) 


where xp has the following values with respect to the different con- 
figurations shown in figure 5: 


In plan form A 


x-l = 1 


( 36 a) 


In plan form B 


In plan form C 


x-> = 1 - £ 

- 1 m 


( 36 b) 


X-] = 1 + •*- 

1 m 


( 36 c) 


In plan form D 


x-, = 1 + 


m 


X 1 2bC 


for 0 g y < m l 2 bC* ’ 1 


for 111(^77 " l) < y ^ ^j~ 


2bC 


( 36 d) 


After the expression for 0, given by the sum of equations (25), (26), 
and ( 27 ), is substituted into equation (35) and the integration is per- 
formed, the results may be reduced to the form 


P = -4pbV 2 k 2 e ;kl)t 


17 (Li + iLg) + + iL 4) 


(37) 


where 


Li = 


x l\/ c2x l 2 


•(4m 2 x 1 - P 2 a 0 ) 


A oy . .1 Cx i 

+ cosh 

C y 


(38a) 
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L 2 = ^x-j 2 - y 2 


T - - 2 |A i) - + 2|A l) 


( 38 b) 


L 3 = \/c2 Xl 2 - y 2. 


^ ^6p 2 M 2 X, 2 - pVi - 31^0! + 2M 2 3 2 x 1 j + 


_4£ 


3(3 2c 2 


P 2 a-i_ - 3P 2 M 2 C 2 A 3 _ 3M 4 c 2 a2 _ 2M 2 X,^j 


* - 2 *0 L 1 


( 38 c) 


34 = \ Jc 2 * i 2 - y 2 


|~(2p 2 A! + P 2 Aq . 4M 2 Xi) + ^ ^ l3 (2pV 2 


P^k 


M‘ 


■p 2 a 1 + 4m^ct 3 - 4p 2 M 2 a 5 ) - ~p^T (k$\: 2 Ao + 2p 2 \ 2 - M 2 a } 


4M 2 p2c 2 a2 _ 8M^C 2 a4 - 8p2c 2 M 2 CTgj 


- cosh 


_1 2 * 2 - 


Acr 


kC 


y^(4p2c 2 A 3 - 2p 2 X 2 + U 2 a ± + 4M 2 p 2 C 2 a 2 ) 


M 2 : 

?C3 


- 2XQL2 (38d) 


In a similar manner, the moment (positive leading edge up) about 
the axis x = xq is 

^1 


/ •"•X 

(x - xq ) Ap dx 

y/c 


= -4pV 2 k 2 b 2 


- (M x - iMg) + o 0 (M 3 + iM U ) 


(39) 
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where 


M- 


1 = /c 2 x 1 2 - y 2 


4x-l 2 


3P‘ 


(hu\ - 6%) + + e 2 A 0 ) 


- 2x 0 L 


Mg = \fc^x^ - y 2 


0 ^^l 3 / o 2 \ 

~ X 1 + — —(& \ - SM 2 ^) 


M^kx- ^ 

p\? 


+ M 2 cr 1 + 4M 2 p 2 C 2 a 2 j 


V 2 


y^P 2 ^ + M 2 a^ + 4M 2 C 2 CT2j 


Ck 


Cx x 

cosh - 2 xqL 2 

y 


Mo = \Jc 2 x 1 2 - y 2 


(p^ a 1 “ 2M 2 P 2 X, 1 - 6M 2 P 2 X 2 + 3M^(JjJ 


-^-g^pV-L - 4M 2 P^C 2 A 3 - 2p 2 M 2 X 1 + 2M 2 p 2 \ 2 - M^ct-l - 4mV 2 C 2 ct 2 
3 C 


A o y2 y^ 


+ -J-j (pV-L - 4M 2 p 1+ C 2 A 3 - 2M 2 p 2 X 1 + 2M 2 P 2 X 2 - M^ - 


Ck 2 pV 3 


4mV 2 C 2 a 2 ) 


C x 

1 -y- - 2x o( M l + L 3 + 2x 0 L l) 


cosh - 


(40a) 


(40b) 


) * 


(40c) 
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= \/c 2 x 1 2 - y 


4x x 2 


3P 2 k 


(«% + ® 2a i - 1,M %) - ^i^( p2A o - » 2a i 


aA-J + 16 m ky -1 _ p2 M 2 0l + kH h O-_ - lm% 2 0 j'j - 


l6M^kxp^y^ 

15P 6 c 2 


^5P 6 C 2 A 3 + 2pV 2 - M 2 P 2 o 1 + 5M 2 pV) 2 a 2 - 


10M^P 2 C 2 a^ + M 2 P 2 a<- - 10M 2 ^C 2 a^j + (p^C 2 A^ - 


4p\ 2 + 2M 2 p 2 a-L + 5M 2 p^C 2 a 2 + 2M^a 3 + 5M 1 j ‘P 2 C ^ 


2M 2 P 2 ct^ + 5M 2 P^C 2 agj - 2x0^^ + + 2xQL^j 


(40d) 


In equation (37) > for example, the quantity (L-^ + il^) is the lift-force 

coefficient associated with vertical motion of the wing. The real part 
is in phase with the vertical position of the wing and the complex 
part L 2 is 90° out of phase with this position. Similar definitions 
apply to the lift coefficient (L^ + iL^) associated with pitching motion 
and to the moment coefficients (M-, + iM^) and (M^ + iM^) . The complex 

or out-of -phase terms determine the aerodynamic damping associated with 
different wing motions. 


Although the expressions for the components of lift-force and moment 
coefficients in equations ( 38 ) and (40), respectively, are lengthy, they 
may he quite easily evaluated with the aid of the graphs in figures 3 
and 4. 
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DISCUSSION 
Sample Calculations 


In order to give some indication as to the general nature of the 
spanwise distribution of the different components of lift and moment 
coefficients, equations ( 38 ) and (40) have been evaluated at different 
spanwise positions y for plan form A (fig. 5(a)) for the following set 

of conditions: PC = 0.5> C = 1.0, Xq = 0.6, M = ^5/4, and k = 0.1. 

These sample results are plotted as functions of spanwise position in 
figure 6. The spanwise variations of the different components of lift 
force are shown in figure 6(a) and the corresponding variations of 
moment coefficients in figure 6(b). 

In figure 6 note that, for the particular set of conditions for this 
example, the maximum values of the components of moment coefficients M^_, 

M 2 , and M 3 are positive and act near the tips of the wing, whereas the 
maximum value of the component M^ is negative but also acts near the 

tips. It may also be noted that the integrated (in spanwise direction) 
values of the components of moment coefficient, or components of total 
moment coefficient, would in each case have the same sign as the maximum 
value of the corresponding component of section moment coefficient. This 
result is not necessarily true in general, because changing some of the 
parameters involved in the evaluation of the spanwise distribution of 
some components of both force and moment coefficients may change the 
distributions significantly from those shown in figure 6. 

The fact that the total component of moment coefficient M^ is 

negative in the example just discussed shows that, for the conditions of 
the example, this term would not contribute to the aerodynamic damping 
but, on the contrary, would act as a source of energy for the oscillating 
system. This circumstance is significant since it leads to the possi- 
bility of the single-degree-of-freedom torsional instability discussed 
in the following paragraphs. 


Undamped Torsional Oscillations 

The wing plan forms discussed herein, like two-dimensional and 
rectangular wings, exhibit the possibility of undamped torsional oscil- 
lations for certain ranges of Mach number M and location of axis of 
rotation xq. This fact is borne out, as indicated in the preceding 
paragraph, by considering the integrated (spanwise) value of the com- 
ponent of damping moment M^ associated with pitching or torsional 
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motions. The main results of this phenomenon can he obtained by con- 
sidering very slow oscillations so that only terms in equation (40d) 
for involving the reduced frequency k to the order l/k need be 

retained. In this case, 




8(2M 2 - 1) 
+ 3P 2 C 2 


*ay 2 


+ 



(2M2 - l)A 1 )x 1 + PSxqAq * 


(kl) 


For the triangular plan form (fig. 5(a)) the integrated value or total 
component of moment coefficient is 


77 110 

M 4 = _ 
2(3^k 


-( 2M 2 + 1)A 0 + 3(2M 2 - 1)A 2 + ^XqJAq - (2M 2 - 1)aJ 


4xo 2 (M 2 - 1)Aq 


(^ 2 ) 


In general, the condition of torsional stability or instability 
depends on the sign and magnitude of in equation (42). Positive 

values of indicate stable conditions and negative values indicate 

the possibility of torsional instability. Between the stable and unstable 
conditions, that is when M^_ vanishes, a borderline state of unstable 

equilibrium separating damped and undamped torsional oscillations exists. 

The ranges of valuesjcf Mach number M and location of axis of 
rotation xq for which Mi^ vanishes for some selected values 

of C = tan € are shown plotted in figure 7* The regions inside the 
curve in this figure indicate instability. The dashed curve, on which 
some of the solid curves terminate, represents the locus of values of M 
and Xq for which vanishes for the whole class of triangular wings 

with supersonic leading edges, that is, for triangular wings where 

PC ^ 1. It will be noted that M4 vanishes for values of Xq ahead of 

the root - 2 -- chord position. It will also be noted that, as the vertex 
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angle e = tan - ^ C decreases to 30°, the range of values of Mach number 
for which vanishes decreases sharply. 

In conclusion, investigation of equation (42) shows that, for a 
given value of the reduced frequency k, Mach number M, and location of 
the axis of rotation xq, the magnitude of the damping coefficient M^ 

decreases as C decreases and, consequently, torsional instability is 
less likely to occur with slender triangles than with wider triangles. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., June 19, 1951 
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APPENDIX 


METHOD FOR DETERMINING DISTRIBUTION FUNCTIONS 


In this appendix a method of obtaining the distribution functions 
is developed in detail. Expressions for the distribution functions 
required to derive the velocity potential to the third power of the 
frequency of oscillations and a set of functions useful in writing the 
expression for the potential in concise form are given. 

As indicated in references 8 and 9 > where the triangular wing is 
treated for constant angle of attack and for constant rolling and pitching 
motions, a convenient form of the velocity potential corresponding to a 
distribution of vertical velocity proportional to x n is 




r 

F n (c) dcr / 


l n+2 as 


\/(x - |) 2 - P 2 (y - ct£) 2 - 3 2 z 2 


(Al) 


where a is the slope of a ray passing through the vertex of the wing 
(n = cr|), is the least value of £ that causes the denominator R 

in the ^integrand to vanish, and F n (cr) = j D n (£,q) is the distribution 

function that is to be determined so that 


( w n) 


z=0 





(A2) 


In equation (Al), the integrand is noted to be singular at the 
limit | = A form of the integrand which avoids this difficulty is 

obtained by making the following change of variables: 


q 



(1 - p 2 ad) 


N = 


_gx 
1 - P : 


- 2 f^T^T7^) 


(A3) 


0 = * 
x 


T = cosh' 


■i q - j 

N 


where t is the new variable of integration. With these substitutions equation (Al) becomes 


0 -ii. 


F n (^) 


cosh - -*- ^ 
^ N 


51 ^-C \jl - p 2 a 2 ^0 


(q - N cosh T) n+ 2 <Jt 


U4) 


,d0 n \ 

and the corresponding expression for w n = - — , which is an integral equation for FJa), 

\° z / z=0 
is 


1 lim 

w n " n z — * 0 ^ 


-p 2 (n + 2) 


' C F n (o) 


-C 


(- 


3 2 a 2 )3/ 2 Jq 


.cosh"- 1 - § . - 

N n ^ - cosh t j cosh t dT + 


P^z 2 (n + 2) ^ 
'■'-C 


C F n (g) do P 

(l - P 2o 2) 5 / 2 J 0 


cosh -1 ^ 

" » n - 2 i 


cosh t ^ ^ + n cosh t ) cosh t dT 


n 


= x 


(A5) 


In this equation the value of w n obtained by performing the indicated integrations and then 
going to the limit z = 0 is the same as the value that would be obtained by first going to 
the limit and then performing the integration but neglecting singularities, pointed out sub- 
sequently, that arise when the value of a approaches the value of 6. Making use of this 
fact reduces the calculations involved and simplifies the integral equation for F n . Thus 

passing to the limit gives for equation (A5) 
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W, 


n 


= g 2 (n + 2) f C F n (q) dq f C0Sh Nq n /j^ 
" J-C (l - pVp/2 j 0 “° L 


cosh cosh t dT 


n 


= x 


(A6) 


where Nq is the value of N at z = 0. The requirement that the normal 

velocity be proportional to x n and independent of y implies that all 
derivatives of equation (a 6) with respect to y or 9 must vanish, or 
that F n (cr) must be so determined that the final value of the integral 
in this equation be independent of 9. The requirement that all deriv- 
atives of equation (A6) with respect to 6 vanish leads to other equa- 
tions for F n and, after n + 1 such differentiations, the equations 

acquire forms for which solutions are known. The value of F n can then 
be determined from these known solutions by evaluating each of the 
d K w n 

derivatives — — at any arbitrary value of 6 in the range -C <9 < C 
d6* 

as follows: 

The Kth derivative of w n with respect to 9 (equation (A6)) gives 
the following integral equation for F n : 


a K w n 

de K 


P 2 (n + 2) 

7t 



Fn(cr) dcr pcosh" 1 ^ / px x K 

(l - (3^2 p/ 2 \1 - 3 2 ct 2 / 


No cosh r) 11 " 1 '-'- i|)(t) dT = 0 

U ' N 0 


(A7) 


\ 
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where 


3>(t) = 


When K 
hinomial 


(n + 1)1 / 3 

(n + 1 - K) ! \N 0 


o + cosh T 
cosh t - q 


k K 


K(n + 1): 

(n + 2 - K) ! 


Per + cosh t 


\N 0 cosh t 


Nr 


+ 


k(k - l)(n + 1) : 

( $(J + cosh T Y 

(n + 3 - K): 

\N 0 cosh T - q ) 


. + (-i) k k:n 0 


-K 


: n + 1, the expression for <D(t) may be recognized as being a 
expression, namely 


$(t) = (n + 1) I 


Pa + cosh t 
Nq cosh t - q 



(n + 1) :(paN 0 + q) n+1 
(~No) n+ l(q - Nq cosh t) 11+ '*' 


(n + l) ! x n+1 

(-No) n+1 (q - Nq cosh T) n+ 1 


(A8) 


Thus the integral equation corresponding to the n + 1 derivative of w n is 


ro 

oo 


d n+1 vn (-l) n+1 (n + 2)10 n+3 x 2n+2 T C Fn( g ) 


ae 


n+l 


J -C 


NQ n+2 (l - P 2 ^ 2 ) 2 


4 u-l x 

J 0 sinh — 

r », 

0 


N 0 \/l-p2a2 


cosh r dr 


(-l) n+1 (n + 2):x n ^/l - p 2 ^ P C F n( q) d0 = 0 

« J-C (0 - a) n+3 


(A9) 


Further differentiation of w n leads to other equations involving integrals similar to that in 
equation (A9) and, as will he seen subsequently, is not necessary for the determination of F n . 

The singularity at a = 6 in equation (A9) is a result of going to the limit z = 0 after 
equation (A5) and, as previously implied, is to be ignored. 

Consider the following equation similar to equation (A9) : 



f m ( g ) dg 

(0 - d ) n+3 


(A10) 


It is known by analogy with problems in incompressible flow and may be shown by direct substi- 
tution and reduction that this equation is satisfied for any value of n (n = 0, 1, 2, . . . ) 
by the function 


f 0 (o) = \/c 2 - cr 2 


(All) 
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This function also satisfies equation (A9) for any value of n and 

satisfies equation (A7) for all derivatives — — beyond the nth 

d0 K 

derivative. A more general form of solution to equation (A10) may be 
shown to be 


f m = o^C 2 - a 2 


(A12) 


where m is an integer and m ^ n. The validity of this solution 
follows, since 


- (0 




m!0 


m-r 


(m - r) !r ! 


(0 - a) r 


(A13) 


and 


o^C 2 - a 2 = mIe m-r / C 2 - a 2 

(0 - cx) n+3 (m - r) Ir ! (0 - a) n " r+3 


(Alh) 


Each term of the summation in equation (Al4) is found to have the form 
of the integrand in equation (A10), from which the function (A12) is 
concluded to satisfy equation (A10) and equation (A7) for all values of 
K > n. 

From the foregoing discussion and consideration the distribution 
function F n (cr) may be uniquely determined in terms of expression (A12) 

as follows: 

Consider the expression 


F n (q) - A 0 f 0 + A^fj + . . . + A n f n (A15) 

where the coefficients Aq, . . . A n are constants that are to 

be determined. 

Each term on the right of equation (A15) is noted to satisfy equa- 
tion (A7) for all derivatives beyond the nth and a total of 
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n + 1 parameters is to be determined. If expression (A15) is intro- 
duced into equations (a 6) and (Aj) and the indicated integrations in 
equation (A 6) and in equation (A7) for K = 1, 2, . . . n are performed, 
n + 1 linear equations are obtained in Aq, A lf . . . A n from which 

these constants may be determined. 

The integrations with respect to a in equations (a 6) and (A7) are 
in general difficult and tedious to perform; however, as previously 
pointed out it is only necessary to perform the integrations for some 
particular value of 0 in the range -C <0 < C. The integrals have 
their simplest form when 0=0 and the integrations can be made for 
this value of 0 by reductions and use of formulas in reference 11. 

The distribution functions F 0 ( cr) , F x (a), F 2 (a), F^cr), Fq( a) , 
and F x (0) and functions D n (|,T|) calculated by the foregoing method are 
as follows: 














where 




P^C^F' + (1 - 3 2 C 2 )E 
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a 2 = 


a 3 = 


A k = 



A 6 = 



A 8 = 
*9 = 


2d 3 


2dxd 3 

+ P 2 C 2 d 2 2 


2d2 

2d!d 3 

+ P 2 C 2 d 2 2 


6d6 

6di(.d6 

+ P 2 C 2 d 5 2 


6d 5 

£'h d 6 

+ P 2 c 2 d 5 2 

2d 3 a 7 

- P 2 C 2 d 2 dQ 

2d x d 3 

+ P 2 C 2 d 2 2 

2 ( d l d 8 

+ 

2d x d 3 

+ p 2 c 2 d 2 2 

6d 6 d 9 

- P 2 C 2 d 2 d 5 

6d^dg 

+ P 2 C 2 d 5 2 

6( dgd^ 

+ d^d^) 


6d k d 6 + p 2 C 2 d 5 2 


The quantities d n are defined as 

= 5(3 2 C 2 - + (2 - lOp^ 2 + gp^C^E' 

1 2(1 - p ^ 2 ) 2 
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2P 2 C 2 F' - (l + (3 2 C 2 )E' 
(1 - p2c 2 ) 2 


(33 2 C 2 - 5pV^)F' - (6 - lOp^ 2 + 2p^C^)E 
(1 - p^ 2 ) 2 


(2732c 2 - 313^:^ + i2p6c6)F' + (6 - 553^2 + - 24 p 6 c 6 )E' 

6(1 - p 2 c 2 )3 

(93 2 C 2 - p^^F' - (3 + 7p2c 2 - 2p^c^)E' 

(1 - 3202)3 


( 3p2c2 - 93^ - 2p 6 C^)F’ - (6 - 193 2c 2 + - ^p^C 6 )E' 

(1 - p2c 2 )3 

* 


3 2 C 2 (E' - F') 

2(1 - p2c 2 ) 


p2c2p-' _ E' 

1 - p2c 2 


p^g^P^ 2 - 3)F' + (U - 2p 2 C 2 )E'] 

6(1 - p2c 2 ) 2 
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Combinations of these functions useful for writing the potential 
in concise form are: 


A ° - Al = Hi 
2 1 

H " 

A 1 - A 2 

2 - k 2 

a 3 

2 

A 2 " a 4 
2 3 

X 1 


3 

& 

II 

$ 

1 

£ 

Xq 

2 

& 

1 

>> 

UJ 

II 


H ■ x 6 = H 

x 6 ■ 

$ 

II 

«> 

+ 

£ 

+ 

A 0 A 6 , 

2 - *9 


- Mo 




a 1 a 9 

2 = M.2 


X9 - Xu = X13 



H . 

M 2 

X 10 


M‘ 


Vf 




= CT 2 


a 3 


= CT| 


*•11 

M 2 


= ap 


+ H2 = 




a 6 


Ho ■ H 2 = x i4 


It will be noted that the quantities Aj_ (i = 0, 1, ... 9) and 
Xj (J - 1, 2, . . . 14) are functions only of the product PC . They 

are shown plotted in figures 3 and 4, respectively. The quantities 
cr k (k = 1, 2, ... 6) are functions of Mach number M and the 

product PC and can be evaluated for particular values of M from the 
plots of A± and \ j . 
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Combinations of these functions useful for writing the potential 
in concise form are: 


2 

A 2 - A^ 


2 A 3 - 


: X 1 

b 

1 

1! 

A, 2 

a 3 

X 10 _ 


2 

M^ - 

x 3 

x 7 

_ hi m 


3 

M 2 

— Xlj. 



x 5 

\Q 

2 

X X 4 

M 2 



_ ^11 = 

*7 

x 6 

M 2 

*6 


II 

col 

£\%: 


* 0*6 , 

— - X, 

i^S-Xu 

a l A 9 

~Y~ = *-12 

X9 - \11 = X13 

*10 " H 2 = ^14 


= a 2 


5 

a 6 


It will be noted that the quantities Aj_ (i = 0, 1, ... 9 ) and 
Xj (J = 1, 2, . . . 14) are functions only of the product 3C . They 

are shown plotted in figures 3 and 4, respectively. The quantities 
a k (k=l, 2 , . . . 6 ) are functions of Mach number M and the 

product PC and can be evaluated for particular values of M from the 
plots of A^ and Xj. 
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Figure 1 



(a) Plan form (xy-plane). 



(b) Section y = y^ (xz-plane). 

- Sketch illustrating coordinate system and the two degrees of 
freedom a and h. 



Figure 2.- Sketch illustrating region of integration for the velocity 

potential. 
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(a) Plan form (xy-plane). 



(b) Section y = (xz-plane). 

Figure 1.- Sketch illustrating coordinate system and the two degrees of 

freedom a and h. 
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Figure 2.- Sketch illustrating region of integration for the velocity 

potential. 
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Figure' 3-- 


Variation of the 


quantities 


Ai 


as functions of 


pc. 
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(a) j = 1 , 2, and 3. 

Figure k.- Variation of the quantities Xj as functions of PC. 
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(b) j = 4, 5, 6, 7, 8, and 9. 
Figure 4.- 


Continued 
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(c) j = 10, 11, 12, 13 , and l4. 


Figure 4.- Concluded 


05 



(b) Plan form B. (d) Plan form D. 

Figure 5.- Sketches illustrating different plan forms for which the force 
equations ( 38 ) and moment equations (40) apply. 
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(a) Lift force. 


Figure 6.- Sparrwise distribution of components of lift force and mome nt 
coefficients for Xq = 0.6. pC = 0.5; C = 1.0; k = 0.1; M = ^5/4. 




(t>) Moment. 
Figure 6.- Concluded. 
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Figure 7.- Curves showing ranges of Mach number M and axis of 

rotation Xq for which the aerodynamic torsional damping moment 
vanishes for some selected values of angle e. 
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